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Abstract Cells constantly need to monitor the state of the environment to 
detect changes and timely respond. The detection of concentration changes of 
a ligand by a set of receptors can be cast as a problem of hypothesis testing, 
and the cell viewed as a Neyman-Pearson detector. Within this framework, 
we investigate the role of receptor cooperativity in improving the cell’s abil¬ 
ity to detect changes. We find that cooperativity decreases the probability of 
missing an occurred change. This becomes especially beneficial when difficult 
detections have to be made. Concerning the influence of cooperativity on how 
fast a desired detection power is achieved, we find in general that there is an 
optimal value at finite levels of cooperation, even though easy discrimination 
tasks can be performed more rapidly by noncooperative receptors. 


Keywords Sensing • Cooperativity • Hypothesis Testing • Stochastic 
Processes 


1 Introduction 

The ability to acquire, process and stock information is crucial for living 
beings. Already at the level of a single cell the external environment is chem¬ 
ically sensed and the acquired information is processed to prepare a suitable 
response. All the steps involved in these processes are permeated with noise: 
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the arrival of the external stimulus through diffusion is a stochastic process 
and so is the interaction with the receptors. The successive biochemical steps 
for the signal transduction usually involve a limited amount of proteins so 
that fluctuations play a prominent role. In other words, cells need to extract 
useful information from noisy inputs through noisy mechanisms, a seemingly 
daunting task. Despite these hurdles, cells are able to sense the environment 
and make decisions with remarkable precision. For a proper investigation of 
sensing it is fundamental to precisely characterize it and to assess its per¬ 
formance. The first step towards a quantitative understanding of this issue 
was taken by Berg and Purcell in their seminal paper in 1977 [6] addressing 
the question of how precise can a cell be in determining the concentration 
of an external ligand. Several studies have followed their line of reasoning 
and considered signal to noise ratios as a measure of the quality of sensing 
(see e.g. p!|l^l8ll9 l[l7[ll8ll23[l24[l34l[35ll42[l46ll48ll49ll5Q[l51j b A complementary 
approach within the framework of information theory has been taken, in¬ 
vestigating the ability of signaling networks to acquire information about 
the environment and transmit it downstream (see Ref. [7] for a general dis¬ 
cussion, [nm for recent reviews and [Tr)l[CT[THll29UdlllHSllHT?lEHllB2U54ll55[ 
for a representative though necessarily incomplete list of re¬ 
cent contributions). Furthermore, the influence of receptor cooperativity has 
been studied (see e.g. p1[^ [^[^l49[l5Qj ). Quantifying the ability to sense, 
transmit signals and respond to the external environment enables a quanti- 

[4^l49[l57j . This contributes to the understanding of the trade-offs which, un¬ 
der evolutionary pressure, may have shaped the signaling and transcriptional 
strategies which are presently observed. 


Here we focus on a different task: rather than precisely inferring the value 
of an external concentration, cells have to efficiently detect changes from a 
reference level. Such systems include for instance the ones involved in en¬ 
forcing homeostasis which must monitor deviations from the physiological 
conditions, the early stages of the immune response with the detection of 
an antigen and more generally signaling pathways downstream of receptors 
that undergo adaptation to the external stimuli. The problem of detecting 
a change can be cast in terms of hypothesis testing. Namely, one tests the 
hypothesis that a change has occurred vs the null one. This problem has 
been thoroughly studied in Ref. m, where it has been shown how the oc¬ 
cupation history of a single receptor can be used to perform a sequential 
probability ratio test. In such dynamic formulation they derived how quickly 
a decision between two hypotheses can be made given an error threshold. To 
carry out a statistical test exploiting the history of a receptor, cells need to 
devise molecular strategies to record it and analyze it. Refs. ESlllT] suggested 
some possible mechanisms to encode statistical analysis in the level of some 
readout molecules which can be used to make decisions. At variance with this 
approach, which focuses on sequential testing and considers a single receptor, 
here we study the case of parallel testing through the instantaneous occupa¬ 
tion state of a large pool of independent receptors. This static approach no 
longer requires the presence of additional molecular layers as the receptors 
state at a given time can be used directly as a readout. As a shortcoming. 
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restricting to the instantaneous receptors occupation provides a less powerful 
test than the one exploiting the full receptor history. 

We then view the cell as a Neyman-Pearson detector which compares 
the likelihood of the observed receptor occupancy distribution under the 
two hypotheses: the environment has not changed vs a change occurred. 
Within this framework we quantify the detection performance in terms of 
the probability of missing the detection of a change. 

The specific question that we address is whether cooperativity between 
the binding sites of a given receptor can be beneficial for effectively detecting 
changes. We find that cooperativity indeed increases the detection sensitiv¬ 
ity making it the preferable mechanism when difficult detections have to be 
performed. When considering the time needed to achieve a certain statistical 
power we unveil a trade-off which results in an optimal finite level of coopera¬ 
tivity which depends on the required sensitivity and on the number of binding 
sites present in each receptor. In summary, we find that easy detections are 
often achieved more rapidly by receptors with low levels of cooperativity (or 
even noncooperative binding sites) whereas for difficult discrimination tasks 
a high cooperativity is mandatory. 


2 Neyman-Pearson hypothesis testing for concentration 
discrimination 

To illustrate the main ideas of our approach let us start with a simple ex¬ 
ample. Consider the case in which a cell has to determine whether the con¬ 
centration of a given ligand has changed or not by means of the occupation 
state of its receptors. Such problem can be addressed in terms of testing 
the null hypothesis that no change occurred vs the alternative one that the 
concentration has changed. For the sake of simplicity let us restrict to the 
specific case in which the concentration has been constant for a long time 
and the change is an instantaneous switch from a value c to c' taking place 
at time t = 0. If a cell is equipped with 1 independent receptors each 

with L different binding sites (see figure Q, the state of the system at a 
given time is specified by ^ ^ 2 , • • • where each = 0,1,... L denotes 

the number of occupied binding sites in receptor n. Since the receptors are 
independent we have that .C is a collection of N independent identically dis¬ 
tributed variables drawn from the probability of the occupation number of a 
single receptor pt(0- Inferring if the concentration of the external ligand has 
changed corresponds to choosing between two hypothesis 


Ho : pt{l) = Po = P%q{l) (1) 

H^ : pS)=Pi=pi{l) 

where Peg(/) is the equilibrium probability of a receptor having I occupied 

binding sites when the concentration is c and pf (/) is the time dependent one 
for a system that started in equilibrium with c but for which the concentration 
switched to c' at t = 0. Neyman-Pearson lemma ensures that for a given 
significance a (probability of mistakenly rejecting the null hypothesis when 
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N=i4 m 
L = 4 Y • 



Fig. 1 Schematic view of a cell and its sensing components. The cartoon depicts 
a case with = 14 independent receptors each equipped with L = 4 binding sites. 


it is true, type I error) the most powerful test (i. e. minimizing the probability 
P of discarding the alternative hypothesis when it is the correct one, type II 
error) is the likelihood ratio [ 12 ] . Such test implies that for a single receptor 
the null hypothesis is rejected when I is such that 

Ail) = < V (2) 

Pt (0 

where 77 determines the significance 

a = ProhiA < r^\Ho) = Y.p%il)9 {ppi{1)-p%,il)) (3) 

I 

where 0 is the Heaviside step function (i.e. 0{x) = 1 if x > 0 and 0{x) = 0 
if X < 0). According to Stein’s lemma, when the number of independent 
receptors is large {N ^ 1) the probability of not detecting the change when 
it has occurred decreases exponentially as P ^ ^-NDKLipoWpi) where the 
Kullback-Leibler divergence is defined as 

OkUpoWpi) = ( 4 ) 


which for this specific example reads: 

DRUplgM') = • 

I Pt 


( 5 ) 
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Fig. 2 Time evolution of the Kullback Leibler divergence Di 
the associated miss probability /3 for iV = 1000 receptors. The level of the divergence 
sets the detection sensitivity. The shaded blue area refers to levels of sensitivity 
which can be reached by the system. As discussed in the text, the concentration 
change cannot be immediately detected as the initial divergence is zero. As time 
goes by, the discrimination power increases until an equilibrium value is reached. 
If one is interested in a specific sensitivity (thin black line in the plot set at a 
miss probability of 5% corresponding to Dkl = 0.003) there is an associated time 
needed to cross the threshold t*. The plot is a numerical solution for the Pauling 
model described in section for a receptor with L = 3 binding sites cooperating 
with a coupling of J = 1. The concentration change to detect is c — 0.9 * Ci/ 2 . 

Time is expressed in units of (fc^^c'/ci/a) i.e. about the time needed for a 
(un-)binding event in the noncooperative case. 


The divergence, and consequently the test sensitivity, changes with time (see 
fig-d- At the very beginning the two distributions are equal and discrim¬ 
ination is therefore impossible. For short times the likelihoods of the two 
hypotheses are still very similar and the Kullback-Leibler divergence grows 
quadratically in time (due to probability conservation there is no linear con¬ 
tribution). For long times the probability of the alternate hypothesis ap¬ 
proaches the equilibrium value limt^oo Pt = Peq divergence saturates 

at the value DKL{p%q\\p%q)‘ The existence of a stationary distribution en¬ 
sures that the Kullback-Leibler divergence grows monotonically in time. This 
means that the asymptotic value DKL{plq\\plq) sets the maximum detection 
precision attainable for the given concentration change. For intermediate 
times, the divergence changes concavity (at least once) from positive to neg¬ 
ative. 
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3 Pauling model of cooperative receptors 

Let us now take a closer look at the probability distribution of the occupation 
level of an individual receptor and study its dependence on the external 
concentration. We consider each receptor to consist of L binding sites and its 
read-out to be given by its number of occupied binding sites I = 0, 1 ... L. 
To investigate the role of cooperativity we choose the Pauling model in which 
each binding site of the receptor interacts with all the other ones [38] . A given 
binding site i = 1,... L, can be occupied (cr^ = 0, 1) with a probability that 
depends on the external concentration of the ligand, on its binding energy 
and, through cooperativity, on the number of other sites of the receptor 
which are bound. The system can be described in terms of an Ising model 
with Hamiltonian: 

^ ^ j 

H =-h'^Gi- (6) 

i=l i^j 

where, for the sake of simplicity, we have encoded in h both the contribution 
of the binding energy and the one of the chemical potential which is affected 
by the ligand concentration: 


h = \og-^ (1) 

where is the dissociation constant for the noncooperative receptor and 
we have set ksT to unity for the rest of the paper. We refer to [38] for a 
detailed connection between the statistical mechanical and the chemical de¬ 
scription of the systems. Notice that, in this system, cooperativity is encoded 
in the fact that the probability of binding increases with the number of bound 
sites. Such cooperativity grows with the number of interacting binding sites 
and with the coupling parameter J. Setting J = 0 corresponds to consider¬ 
ing the noncooperative case of independent binding sites. The energy of a 
configuration depends only on the occupation number I = giving 

H{1) = -hi - - 1) (8) 


and the number of configurations with the same occupation number is simply 
given by the binomial coefficient (^). The equilibrium occupation probability 
then reads 


(f) exp - 1)] 


( 9 ) 


and the average number of bound sites can be directly computed as 


L 


{l) = Y,lPeS)- 

1=0 


( 10 ) 


The cooperativity of the system allows to have sharp changes in the occupa¬ 
tion probability as the external concentration is varied. The system is then 
said to be ultrasensitive around the value of the external field for which half 




7 


of the binding sites are bound on average hi /25 which we refer to as the tran¬ 
sition point. Such value depends on the number of total binding sites L and 
the coupling J 


K/2 = ~{l-i) 


( 11 ) 


and ensures that Peq{l) = Peq{L — 1). We note that for highly cooperative 
systems the equilibrium distribution at h = hi /2 is concentrated at the ex¬ 
tremes I = 0 and I = L. The transition point occurs at different values of 
concentration as the number of binding sites or the coupling are changed 




( 12 ) 


To proceed to a meaningful comparison between different cooperativities 
we need to consider systems that have to detect changes of the same rel¬ 
ative amplitude c/ci/ 2 . In general, cooperative receptors with different num¬ 
bers of binding sites or cooperative strengths can set their transition point 
at the same concentration by tuning their dissociation constant. This is a 
simple mechanism that implements sensorial adaptation. Indeed, by letting 
J) = the transition point is reached for the same concen¬ 
tration Ci /2 = This means that in order to maintain ultra-sensitivity 

around a given concentration the receptor has to compensate for a higher co- 
operativity by evolving towards a higher dissociation constant (lower affinity 
with ligand). 


3.1 Cooperativity and discrimination 

As discussed above, around the transition point, the distribution of the occu¬ 
pation number of the receptor exhibits a stronger dependence on concentra¬ 
tion as cooperativity is increased. Such sharper dependence makes it easier 
to tell if the concentration has changed from the reference value (set at the 
transition point). Intuitively, then, cooperativity improves the discrimination 
ability of the receptor. To check this intuition let us consider the behavior 
of the Kullback-Leibler distance between the equilibrium probability at the 
transition point and the one of a different concentration and its dependence 
on the degree of cooperativity (see figure |^. We start by considering the 
asymptotic value of the divergence which corresponds to the largest possible 
discrimination. Due to Stein’s lemma, the Kullback-Leibler divergence de¬ 
scribes how the probability of missing a change in concentration decreases 
with the number of independent measurements N in the limit N ^ 00 . For 
instance to achieve a miss probability P ^ 5% we need a divergence between 
the two distributions of 



(13) 


For a given number of independent receptors N the Kullback-Leibler diver¬ 
gence then sets the miss probability. 





1 2 3 4 5 6 7 8 


L 

Fig. 3 Long-time asymptotic Kullback Leibler divergence and miss probability 
for iV = 100 receptors as a function of binding sites numbers for c = 0.9 * Ci/ 2 . 
The full purple line is the divergence in the limit of infinite coupling (J oo) 
discussed in eg. (|17| ) and the full blue one is the noncooperative case (J = 0) 
discussed in eq. The shaded blue area depicts the values of precision reachable 
with noncooperative binding sites, the purple shaded one the values achievable by 
cooperative receptors. The white area above the purple line refers to precisions 
which cannot be reached for the given number of binding sites. The cyan points are 
numerically evaluated results for J = 3 already approaching the infinite coupling 
limit. The red points refer to the case with J — 1 and the green ones to J = 0.5. 


Let us start by considering the noncooperative case (J = 0). The system 
then simplifies greatly and is just a collection of L independent two states 
processes with binomial distribution 


pT,{i) = 



(14) 


The Kullback-Leibler divergence is then analytically accessible and consid¬ 
ering the two equilibrium distributions reads: 


T^nc 

^KL 




eq') 


Llog 



where we shall recall that hi /2 = 0. The most relevant feature is that the 
divergence grows linearly with the number of binding sites present in the 
receptor. 

We expect cooperativity to increase the Kullback-Leibler divergence and 
consequently decrease the probability of not detecting a concentration change. 
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In the limit of large coupling J ^ oc we can roughly estimate the occupa¬ 
tion probability by inspecting equation ® around hi /2 = —J(I/ — l)/2 and 
observing that only states with I = 0 ana / = L will have a finite probability: 


giving 


lim p{l = 0) , —r 


-fe) 


(*) 


^L{h'-hi/2) 

lim p(l = L) - TTTJ —^—V — T 


lim DKL{plP\\pin) = log 


J —^OO 


= log 




cosh 


(16) 


(17) 


which for L > 1 is larger than its noncooperative equivalent. Hence, for a 
fixed L, having a cooperative receptor allows to reach values of the Kullback- 
Leibler that cannot be approached by means of noncooperative binding 
sites resulting in higher detection sensitivities. However, the limit in eq. (17) 
cannot be exceeded and this sets an upper bound on the performance of a 
cooperative receptor following Pauling model. 


3.1.1 Small concentration change limit 


The hypotheses we are testing have probability distributions that differ be¬ 
cause of the parameter c and c'. Obviously, when the two concentrations 
are equal the Kullback-Leibler divergence is zero and is at a minimum since 
it is non-negative by definition. Then, for small concentration differences 
Sc/c 1, the divergence is quadratic in the concentration difference. It is 
known that curvature of the divergence in c = c' is equal to the Fisher infor¬ 
mation of the probability distribution with respect to the parameter c which 
is defined as 


3(=) = E 




1=0 


V 


dc 


plS) ■ 


(18) 


Then, by Taylor expansion, the Kullback-Leibler for small concentration 
changes around c reads 


lim 

5c/c^0 


DKLiplMt'") 


\{5cf3{c). 


(19) 


Plugging relation 0 into the equilibrium probability 0 we can compute 
the derivative of the probability with respect to concentration and obtain: 


dp%q{l) 

dc 



( 20 ) 
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which means that the Fisher information is J(c) = ^var{l) and the Kullback- 
Leibler for small concentration changes around the transition point C 1/2 


lim 

dc/c^O 


DKL{Peq^^ 




1 

2 



2 

var{l) . 


( 21 ) 


It is interesting to compare the two extremes of large cooperativity and 
vanishing cooperativity. Taking the large coupling limit J ^ oc of the vari¬ 
ance we see that 


lim lim Oxhiplq^ 

J^oo (5c/c^0 


lb 


Cl/2+<5 c- 

eg 


lim - 

J —^OC 2 



2 

var{l) 


1 

8 





^ (22) 

where two limits can be shown to commute by considering eq. ( p!7| ). For the 
noncooperative case one finds 


lim 

6c/c^0 


'KL 


{p 


Cl/2 

eq 


\\ptr^'^) 


1 / 5c 

8 bi/2 


2 

L . 


(23) 


The strongly cooperative case displays a quadratic dependence on the num¬ 
ber of binding sites, as opposed to the linear one of the noncooperative model. 
The miss probability of the test then decreases exponentially with L‘^N for 
the cooperative case and simply as LN for the noncooperative one. For finite 
couplings, the variance has to be computed directly from the distribution and 
reads: 


var{l) = {l)f Peq{l) (24) 

z=o 

^ Ef =0 (t) exp [hi + - 1)] _ / Ef=o I (t) exp [hi + jl{l - 1)] \ ^ 

Eto (t) exp - 1)] V Eto (t) exp [hi + jiq - 1)] ) 

Such expression, due to the terms exp [^l{l — 1)] does not display a simple 
scaling in terms of L for finite couplings. 


3.1.2 Connection between Hill coefficient and Kullhack-Leihler divergence 


The Hill coefficient is a measure of the cooperativity of binding and can be 
defined (see e.g. Ref. m) as: 


riH = T 


4 d{l) 


L d log c 


4 d{l) 
L~dh 


(25) 


C=Ci/2 C=Ci/2 

Performing the derivative of the mean value with the explicit expressions 
given in eqs. ^ and (10) we see thalQ 


dh 


= var{l) . 


(26) 


^ This is simply the fluctuation dissipation relation between the susceptibility 
and the magnetization variance for ferromagnetic systems 
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Fig. 4 Hill coefficient as a function of binding sites per receptor L for different 
values of the coupling J. For the noncooperative case the Hill coefficient is 1 whereas 
it is L in the limit of strong cooperativity. The parameters and color code are the 
same as for figure 


We can then write 


uh = Yvar{l) 

1j 


(27) 


As we have seen in eq. (19), for small concentration changes, the Kullback- 
Leibler divergence is proportional to the Fisher information which, for the 
system we are considering, is itself proportional to the variance of the dis¬ 
tribution J = var(l)/c^. Hence, for small concentration changes around the 
transition point, we have that 


= \ (*) vard) = | „„ (28) 

linking the Kullback-Leibler divergence and the Hill coefficient. The Hill co¬ 
efficient ranges from 1 for noncooperative binding to a maximum of L ap¬ 
proached in the limit of infinite cooperativity (see fig. 0. This offers some 
insight on the scaling behavior of the Kullback-Leibler divergence for finite 
couplings. Indeed, for small concentration changes, the Kullback-Leibler di¬ 
vergence scales as Luh where uh is the Hill coefficient of the receptor. 


3.2 How long does it take to detect a change? 

From figure we know that at the moment in which the concentration 
changes it is not possible yet to detect the change as the distributions of 
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the two hypotheses are still equal. Detection becomes possible as the like¬ 
lihood of the alternative hypothesis evolves in time. We are interested in 
studying how long it takes before the distributions of the two hypotheses 
are different enough to allow a detection with a given sensitivity. In order to 
do so, let us now consider the kinetic behavior of a receptor whose binding 
sites follow the Pauling model. In general the occupation probability evolves 
according to a master equation 

^ ^ k{l' ^ l)p{l') - k{l ^ l')p{l) (29) 

I 


where, for the system to reach the correct equilibrium, detailed balance must 
be satisfied by the rates: 


k{v^V) p,S) VV^VA 

which for the two allowed transitions implies 

k(l 1) ^ L - I hj^ji ^ {L-l)c . . 

+ / + ! (/ + 1)ci/2(J,L) • ^ ^ 

There is freedom in how to choose the individual rates still satisfying detailed 
balance. We consider the case in which the rates are exponential in the energy 
difference associated with the transition and to make contact with the non- 
cooperative binding we set them to 


k(l I ^1) — {L — l)ku - ( j T\ 


J 

e 2 


0 - 


V) 


(32) 

(33) 


where ku is the unbinding rate for noncooperative binding sites. We can 
then proceed to the numerical solution of the master equation and obtain 
the time evolution of the system. We shall focus on the case in which the 
system starts in equilibrium with the concentration at the transition point 
Ci /2 and the concentration is switched to a different value c'. The most re¬ 
markable feature is that, for fixed d as cooperativity is increased (both 
by having more binding sites or by a higher coupling J) the system slows 
down. The slowing down is more severe for small concentration changes. 
This phenomenon is related to the critical slowing down (discussed also in 
Ref. [48]) and is due to the fact that a system with strong interactions is 
“locked” into a configuration due to cooperativity and reacts more slowly to 
a change in the external field. 


Let us now investigate how this impacts on the time required by the pool 
of receptors to reach a certain detection sensitivity. We need to consider the 
dynamic behavior of the Kullback-Leibler divergence during equilibration to 
the new concentration d and how rapidly it manages to reach the threshold 
associated with the desired sensitivity as sketched in figure As we have 
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b) 


p=5% 



N=1000 
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Fig. 5 Time needed to achieve a desired sensitivity in detecting a concentration 
change to d — 0.9 * Ci /2 for different numbers of binding sites and cooperativity. 

This corresponds to time needed by the Kullback Leibler divergence Dkl{p%^^ \ \pt ) 
to cross a given threshold, a): miss probability of 5% which for N = 1000 recep¬ 
tors corresponds to Dkl — 0.003. For J < 0.81 the receptor with L — 2 binding 
sites cannot reach a divergence of 0.003. For L = 3 the shortest time is obtained 
at finite coupling whereas for L = 5,8 for the noncooperative case. Notice the 
exponential slowing down observed for L — 8. b): miss probability of 0.1% corre¬ 
sponding for N = 1000 receptors to Dkl — 0.0069. Now the required sensitivity is 
unattainable with L — 2 binding sites. It requires a coupling J > 1.37 for L = 3 
and a weak J > 0.1 for L = 5. For the shown number of binding sites the desired 
precision is fastest reached for finite level of cooperativity. Time is expressed in 

units of (^u y^c'/ci/ 2 ^ i.e. about the time needed for a (un-)binding event in 
the noncooperative case. 


discussed in the previous section and shown in figure the divergence be¬ 
tween the equilibrium distribution in Ci /2 and in d is higher for cooperative 
receptors consequently allowing to reach levels of sensitivity which are unfea¬ 
sible for noncooperative binding sites. However, a larger cooperativity slows 
down the dynamics. When considering how fast a certain level is reached we 
are then facing a trade-off as shown in figure 
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To obtain some insight let us study a specific example. Consider for in¬ 
stance the case in which we have N = 1000 receptors, the concentration 
decreases of about 10% from its transition value: c' = 0.9 * Ci /2 and we want 
to achieve a miss probability of /3 5%. This sets the Kullback Leibler dis¬ 

tance between the equilibrium receptor distribution at the transition point 
and the evolving one at Dkl{p%^^ \\Pt) — 0-003. From eq. ( [l^ we know that 
such sensitivity cannot be achieved by a receptor with 2 binding sites unless 
they operate in a cooperative fashion. Then, for small couplings, a receptor 
with 2 binding sites will never reach the threshold. As the coupling is in¬ 
creased the threshold becomes attainable and it progressively takes shorter 
times to approach it until an optimal value is reached. When several binding 
sites can be employed the situation is quite the opposite as the threshold is 
largely exceeded also in the noncooperative case so that the main effect of 
cooperativity is the exponential slowing down of the dynamics. In such case 
the required sensitivity is reached earlier by noncooperative binding sites. 
For intermediate scenarios, as for example with 3 binding sites, the nonco¬ 
operative case is sufficient to cross the threshold but its asymptotic value is 
not much larger. Then, a weak cooperativity is slightly faster as a result of 
the trade-off between a higher equilibrium limit and slower overall dynamics. 
The representative cases of L = 2, 3, 5,8 are plotted in figureLet us study 
how requiring a higher sensitivity with a miss probability as low as /3 0.1% 

impacts the trade-off. We first notice that such low miss probability is not 
achievable if each receptor has only 2 binding sites since N = 1000 noncoop¬ 
erative such receptors can at most reach a miss probability of /d ^ 0.7%. For 
1/ = 3, 5 the required sensitivity can be attained only by cooperative recep¬ 
tors and it is most rapidly reached for intermediate coupling values. Finally, 
employing 8 binding sites cooperativity is no longer necessary to reach the 
threshold. However, a small coupling (J 0.14) allows to cross the desired 
value at slightly shorter times. 


In general, considering how cooperativity affects the time need for a sensi¬ 
tive detection we can identify three main classes which are determined by the 
level of desired sensitivity and the number of available binding sites. Namely, 
what matters is how large the long-time asymptotic discrimination power 
of the noncooperative receptors (eq. [T^ is compared to the one we want 
to obtain. If the required sensitivity is easily achievable by noncooperative 
receptors (threshold much lower than the equilibrium value for the noncoop¬ 
erative receptors), cooperativity is detrimental for a rapid decision (see e.g. 
the case of L = 8 binding sites and [3 = 5% depicted in cyan in figure |^a). 
For more sensitive detections, which are barely attainable by noncoopera¬ 
tive binding sites, a small degree of cooperativity can slightly speed up the 
system as shown for instance for L = 3 and /d = 5% in blue in figure a. For 
hard detections (desired sensitivity higher than the one provided by nonco¬ 
operative receptors) the system must be cooperative. There is a large but 
finite coupling minimizing the time needed to achieved the threshold as re¬ 
ported for 1/ = 3 and L = 5 for a miss probability of /3 = 0.1% in figure |^b 
(respectively blue and green line). 



15 


So far, we have considered the time needed to reach a certain level of 
precision as a function of the coupling intensity J for the case in which the 
reference concentration is set at the transition point c = Ci 12 - This is the value 




Fig. 6 Time needed to reach a certain sensitivity as a function of the coupling 
coefficient J and minus the external field —h — log^. The solid line indicates 
—h — J[L — 1) / 2 which corresponds to the values of the concentration set at the 
transition point that were used to draw figure |^The top figures refer to the case 
with L = 3 and the bottom ones to L = 5. For each value the concentration 
change to be detected is such that c'/c = 0.9 corresponding to 6h = —0.1. Time is 
expressed in units of k~^. 


for which the system is most sensitive to concentration changes. In order to 
extend our reasonings to cases in which the system is not optimally tuned, let 
us consider how the detection time depends on the reference concentration 
(or, equivalently, on the dissociation constant). The results are reported in 
fig-11 For each level of concentration it is possible to determine the range 
of cooperativity that allows for a given sensitivity and the coupling intensity 
J that reaches it in the shortest time. The fastest detections are achieved 
for combinations of concentration and coupling such that the system is close 
to the transition point but higher concentrations are favored. This reflects 
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the fact that the concentration change has to be within the dynamic range 
of the receptor and that higher concentrations imply a faster ligand binding 
consequently accelerating the dynamics. 


3.3 Non-homogeneous coupling strength and ligand affinity across the 
receptors. 


In the previous sections we have considered the case in which each receptor 
has the same coupling strength and ligand affinity. In general, it is possible 
that different receptors have different J and Let us consider the case in 
which the J and of each receptor are drawn from a distribution p( J, Kd). 
If the number of receptors N is large enough so that each value of the pair 
J^Kd is sampled several times we can still exploit Stein’s lemma for each 
receptors subpopulation. The overall miss probability is then an average over 
the different receptors subpopulations reading 

log /3^-N I p{J, Ka) Dkl Ka) \ \pt' {J, K^)) dJ dKa . (34) 

As a specific example, consider the case where the receptors are a mixture 
of two different subpopulations characterized by two different couplings 
and Js with the respective Kd tuned as Kd{J) = so that the 

transition point is for both populations at C 1/2 = 

Let us denote the probability of drawing the strongly cooperative (Jg) 
receptor as q so that we have, on average Nq strongly cooperative receptors. 
The miss probability then obeys 


log^ ■ 


-N 


QDkL {Js)) + (1 - q)DKL {pl,{Jn,)M'(Jn,)) 


(35) 

resulting in a linear combination of the two Kullback-Leibler divergences. By 
varying the parameter q from 0 to 1 the curve smoothly varies between the 
weakly and the strongly cooperative cases. In general, such mixture of recep¬ 
tors displays a faster initial response compared to the strongly cooperative 
case and a higher asymptotic precision with respect to the weakly cooperative 
case (see figl^. This allows for intermediate behaviors which can be tuned 
to satisfy different requirements by changing the sizes of the subpopulations. 


4 Conclusions and discussion 

We have studied the problem of detecting a change in concentration of a 
ligand by means of a pool of receptors. To investigate the detection perfor¬ 
mance we have addressed the problem in terms of hypothesis testing making 
use of the instantaneous occupancy. We have focused on the probability of 
missing the detection of a change that actually occurred and employed it as 
a measure of the quality of the sensing system. With this formalism we have 
considered the influence of cooperativity on the sensing performance. We 
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Fig. 7 In the case of non-homogeneous coupling strength across the receptors the 
miss probability is given by a linear combination of the various couplings. The plot 
refers to the case in which a receptor is cooperative with J — 1 with probability 
q — 1/2 and is noncooperative with probability 1 — q — 1/2 with L = 3 and 


^ Time is expressed in units of (^u y^c'/ci/ 2 ^ 


have identified a twofold effect of cooperativity: a markedly improved sensi¬ 
tivity in the long time limit and a significant slowing down of the receptors 
dynamics. This is consistent with the findings of Ref. [48] in which the sensing 
performance of locally cooperative receptors was considered in terms of the 
precision with which they are able to determine a concentration change by 
time averaging their history. We stress here that both the higher asymptotic 
sensitivity and the slowing down are separately relevant. Most notably, coop¬ 
erativity enables to detect changes with sensitivities that are not attainable 
by noncooperative receptors. When considering the time from the concentra¬ 
tion change needed before the instantaneous receptor occupancy state can 
provide a given detection sensitivity, the two effects of cooperativity trade 
off. In the framework adopted by the authors of [48], the effect of the slow¬ 
down due to cooperativity on the time averaging procedure dominates the 
additional sensitivity and the authors concluded that noncooperative recep¬ 
tors perform better than cooperative ones. On the contrary, we have seen 
that, casting the detection problem in terms of hypothesis testing by means 
of the instantaneous occupancy of a pool of receptors, cooperativity can be 
beneficial. We have found that there is an optimal degree of cooperativity 
which depends on the required sensitivity level and on the number of binding 
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sites present in each receptor. At given number of binding sites we can iden¬ 
tify three general scenarios depending on the required detection power. For 
hard detections the system must be cooperative and the desired sensitivity 
is reached earlier at finite but large levels of cooperativity. For intermediate 
cases the system can achieve the needed statistical power by noncooperative 
binding sites but moderately cooperative receptors are faster at approaching 
it. Easy detections are performed more rapidly by noncooperative receptors. 

It is interesting to consider how the cell’s detection ability scales with the 
numbers of receptors and binding sites. For noncooperative binding sites the 
probability of missing the detection decreases exponentially with the total 
number of binding sites present in all the receptors NL. We have shown 
that for cooperative receptors in the long-time limit the exponential decrease 
is more marked and in the infinite coupling limit for small concentration 
changes depends quadratically on the number of binding sites present in 
each receptor: NL!^. Since it is costly to produce proteins (see Refs.[T5l[T9]) 
sensing mechanisms that improve the performances without increasing the 
number of required components may represent an advantage. The different 
scaling in N and L can then impact on how to invest the limited resources 
favoring either solutions involving few receptors with many binding sites or 
the opposite. The characterisation of the ensuing trade-offs is an issue that 
surely deserves further study in the near future. 
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